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RESUME 


Hydrolab 2018 est un logiciel libre, développé de façon à résoudre les problèmes les plus 
fréquemment rencontrés par les hydrologues. Notre volonté a été tout d'abord de créer un outil très 
simple parfaitement intégré au logiciel EXCEL dont l'utilisation est universelle. 


Les « entrées- sorties » étant des feuilles Excel, il est très facile de les personnaliser avec une 
connaissance minimale de ce tableur. En particulier, il est aisé de modifier la présentation et le 
format des résultats. On peut ajouter des commentaires libre, la langue dans laquelle les résultats 
sont présentés peut être changée en ajoutant une colonne traduisant dans la langue de votre choix, 
la version française. 


Les résultats sous forme graphiques ou de tableau sont immédiatement transférables et 
modifiables sous Word, PowerPoint .... 


Toutes les sources sont fournies sous la forme de macro-commande en basic. Elles peuvent 
éventuellement être modifiées. 


Les principales utilisations concernent : 
e les ajustements statistiques (loi de Gauss, log normale, racine normale, loi de Gumbel, loi 
exponentielle, loi GEV, loi de Poisson) ; 
e la détection d'anomalies dans des séries de données par analyse des résidus de régression ; 
e les régressions multiples et l’analyse en composantes principales ; 
e le comblement de lacunes dans des séries de données. 


Vous trouverez également des utilitaires sur : 
e la reconstitution des crues à partir des précipitations (différentes fonctions de production, 
d’hydrogrammes unitaires et de pluies de projet) ; 
e le laminage des crues à la traversée de barrages (type et nombre d’évacuateur variables) ; 
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AVERTISSEMENTS 


Hydrolab 2018 est un logiciel libre, simple et gratuit. La seule contrepartie demandée est de 
signaler son utilisation et son origine ! 


Hydrolab 2018 peut facilement être modifié mais la contrepartie est que les résultats des 
traitements restent sous la responsabilité des utilisateurs. Nous vous engageons vivement à 
n’utiliser que la version diffusée sur le site de Polytech’ Nice ou sur des sites de confiance. 


Afin de préserver la reconnaissance de l’origine de ce travail, Hydrolab est référencé auprès de 
l'Agence pour la Protection des Programmes 


(CN : 


URL du certificat de référencement : 
http://www.legsalis.net/csi-iddn/certificat.ce1?IDDN.FR.010.0075748.000.R.C.1999.027.20700 


Au fur et à mesure, Hydrolab s'enrichit, se perfectionne et se corrige en fonction des utilisations 
qui en sont faites. Aussi n'hésitez pas à nous faire part des anomalies que vous auriez détectées, 
des améliorations que vous proposeriez et des développements que vous souhaiteriez, à l'adresse 
suivante : laborde@unice.fr 
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INTRODUCTION 


Notre volonté a été tout d'abord de créer un outil très simple parfaitement intégré au logiciel 
EXCEL dont l'utilisation est universelle. Pour se servir correctement d’Hydrolab, il est donc 
préférable d'avoir des bases sur l'utilisation de ce tableur. 


L'objet d’Hydrolab n'est certainement pas de remplacer des logiciels de statistiques beaucoup plus 
complets tels qu'on peut en trouver dans le commerce. Par contre son utilisation est des plus 
simples et des plus rapides. Il présente également une souplesse d’emploi dans la manipulation des 
données d’entrée et dans la présentation des résultats (graphiques et numériques). Enfin il vous est 
possible de modifier certains traitements et/ou de les automatiser. C’est un logiciel ouvert ! 


Avec l’arrivée sur le marché d’'EXCEL 2013, et ses profonds changements par rapport aux versions 
précédentes, il était nécessaire de réaliser une nouvelle version. Ceci a été aussi l’occasion de 
compléter les traitements proposés et de donner plus de souplesse dans l’organisation des données 
en entrées. 

e Dans la plupart des traitements, on sait a priori que les données seront unidimensionnelles. 
Hydrolab 2018 permet dans ce cas de sélectionner les données dans n’importe quelle feuille 
de calcul et de lancer directement le traitement. 

e Il est possible d’importer directement les noms les variables et les observations ou de 
générer automatiquement ces noms, ce qui facilite l’exploitation des résultats. 

e Dans les plages de données, il existe parfois des lacunes, aussi Hydrolab 2018 identifie 
automatiquement ces lacunes (marquées par des blancs ou des caractères 
alphanumériques), les signale et ne fait de traitement que sur les données réellement 
utilisables. 

Dans ce cas de figure il suffit de sélectionner dans votre feuille les données et éventuellement les 
noms des observations et/ou variables à traiter et de lancer le traitement désirer. 


Pour les données multidimensionnelles, Hydrolab 2018 vous impose une mise en forme du tableau 
de données. Comme précédemment il suffit de sélectionner dans votre feuille les données et 
éventuellement les noms des observations et/ou variables à traiter mais vous devrez les coller 
(collage spécial « valeurs » avec ou sans transposition) dans une feuille modèle où vous pourrez 
préciser la nature du traitement désirer et/ou la précision souhaitée. Ici aussi les lacunes seront 
gérées automatiquement. 


Enfin, Hydrolab 2018 comporte des utilitaires spécifiques et détaillés plus loin, sur : 
e la reconstitution des crues à partir des précipitations (différentes fonctions de production, 
d’hydrogrammes unitaires et de pluies de projet) ; 
e le laminage des crues à la traversée de barrages (type et nombre d’évacuateur variables) ; 
e l’évaluation des ETP et ETR selon différentes méthodes classiques. 
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I. 


I.1 


INSTALLATION D’HYDROLAB 2018 


Décompresser le dossier Hydrolab2018 


Généralement vous recevrez Hydrolab2018 sous la forme d’un dossier compressé auto-extractible. 
Si vous avez reçu directement le dossier contenant les documents excel passez à l’étape suivante. 
Suivant les niveaux de protection d’internet vous recevrez le fichier auto-extractible avec ou sans 
l’extension « .exe ». Si elle est absente comme ici : 


gem 


Organiser v Inclure dans la bibliothèque + Partager avec v Nouveau dossier am 09 
ad Musique Vie Nom [ 
Vidéos u 1 
+4 Groupe résidentiel E 1 
B Laborde Z 
Il suffit d’ajouter manuellement l’extension« .exe ». 
Organiser + Es Ouvrir v Imprimer Graver Nouveau dossier Eu 
& os (C) # Nom | 
c HP_RECOVERY (D:) Ji rc 
É} Lecteur DVD RW (E:) m mc 


(Œ) FREECOM HDD (F:) 
cÀ Lecteur BD-ROM (G:) 


B Hydrolab2010 compressfere ] 


EP Fiacceurt vlev 


Un double clic sur Hydrolab2018 compressé.exe vous permet d’installer le logiciel. Attention, 
vous devez impérativement mettre le dossier asua directement sous la racine C:/ 


a Nom 


di Surfer_setup_v10.3.705 
d Thunderbird 

à Other 

de AAA Disque dur JPL 
à e-academy Inc 


e Press Extract button to start extraction. 


+ Use Browse button to select the destination folder 
from the folders tree. It can be also entered 
manually. 


+ the destination folder does not exist. it will be 
created automatically before extraction. 


a p Start MAC.app 
í U rc 
d mac 
R] Hydrolab2010 compressé.exe 
Classeur1.xisx 
v E —— 
15/2015 17:35 Date de création : 13/05/20! 
Mo 


D Dar Can 


Vous devriez retrouver dans votre degi C, un dossier Hydrolab2018 contenant tous les fichiers 


nécessaires. 

Organiser > Inclure dans la bibliothèque ~ Partageravec œ Graver Nouveau dossier v E 

a È OS (C) # Nom Modifié le T 

à re E) Amortissement d'une er. 13/05/2015 15:02 Fe 

D'asuenser Averse -Cruexls 05/20151500 Fe 

B AdwCleaner rs | Gi] Hydrolab 2010.docx 13/05/20151017 D 

Ji Données pluviométriques d'Algérie x S) Hydrolab2010 xs 13/05/2015 09:25 Fe 

-i sa w) Q10 méthodes empiriqu.…  22/ 5 21:05 Fe 

S E) Tab Regmul.xis ‘2015 10:01 Fe 

RES CRU ja Tab Anomaliexls 9/2014 19:51 Fe 

a _— E) Tab Poissonxis 9/2014 17:45 Fe 
Brass | E) Tab ACPaxis 9/2014 16:17 Fe ç 
D Porte () Tab Variogrammel.xis Fe 

FRERES Tab Racinexls Fe 

s Tab Gumbel.xis Fe 

b magai ne Tab Comblement.xis Fe 

1 Program Files (86) E) Tab GES ra 

4 u un Tab Expoxis Fe 

= E) Tab Galton.xls / Fe 

i ro EE) Tab Gaussads f i Fe 

[3] VB5FR.dIl 14/02/2011 1732 Ex 

D temp 8] VB6FR.dll 14/02/20111730 Ex 

Ji Thunderbird a 
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I.2  Autorisez l’utilisation des macros 


Suivant la version d’excel que vous utilisez, referez-vous à l’aide en ligne pour autoriser 
l’utilisation des macros d’Hydrolab2018. A titre d'exemple voici pour Excel 2013. 

1.Cliquez sur l’onglet Fichier. 

2.Cliquez sur Options. 

3.Cliquez sur Centre de gestion de la confidentialité, puis sur Paramètres du Centre de gestion de 
la confidentialité. 

4.Dans le Centre de gestion de la confidentialité, cliquez sur Paramètres des macros. 


5.Choisisez d’activer toutes les macros 
Centre de gestion de la confidentialité ? X 


Editeurs approuvés 

PP Paramètres des macros 
Emplacements approuvés 
O Désactiver toutes les macros sans notification 


Documents approuvés SC à AA 
O Désactiver toutes les macros avec notification 


Catalogues d'applications approuvés © Désactiver toutes les macros à l'exception des macros signées numériquement 
Compléments © Activer toutes les macros (non recommandé ; risque d'exécution de code potentiellement 
dangereux) 


Paramètres ActiveX 


| Paramètres des macros | Paramètres de macros pour les développeurs 


Mode protégé F1 Accès approuvé au modèle d'objet du projet VBA 


Barre des messages 


6.Cliquez sur OK. 


I.3  Chargez le Solveur 


Bien que ce ne soit pas 
indispensable, il vous sera souvent 


Général - fe 
:@ Affichez et gérez les compléments Microsoft Office, 


utile de disposer des fonctionnalités | vvs Congtiments 

du solveur d’Excel. eg Ds {Ou 1e 

Dans le menu « fichier » cliquez sur ||| sn mn Eeo 
Personmakser le Ruban airs Tecir Due Dan ANASINAMALYSSZIUL Complément Excel 


L.ANAYSANATPYBAERXLAM Complément Excel 
o MMKOMKEEROFFRHD.DU  Inipecteur de document | 
snareg\Smart Tag MOFLOUL Action 
Morel IOFFRHD.OLL Inspecteur de document 
..SMeeOrNicel4VOFFRHD.OL inspecteur de document 
.Mer\OMMcelAVOFFRHD.OLL Inspecteur de document | 
.MerlOrncelAIOFFRHD.DLL Inspecteur de document 
Kit d'extension XML 
\edSmart TagINTMAMEOLL Action 
L IAlibramEUROTOOLXLAM Complément Excel 


« option » puis sur « complément » hnair TooiPak - vaa 
et faire «complément excel- 
atteindre » 

Il vous reste alors à cocher la case 


« Complément Solver » 


Barre d'outils Accés rapide Contenu innssble 
| Date pan) 
cn) rues 

En-têtes et pieds de page 

Centre de gestion de la confidentialité Feuilles de caicut cachées 
Lignes et colonnes cachées 
Microsoft Actions Pane 3 
Nom de personne (Contact Outlook] 
Outils pour l'eure 


nn nhnonanon 


cations désactivés 


Macros complémentaires disponibles : 
le == 
[[ | Analysis ToolPak - VBA 

[FJ] Complément Solver | 
F| Outils pour l'euro 


Complément: Complément Sober 

Éditeur 

Compatibilité: Aucune information de compatibilité dsporible 

Emplacement:  CProgrem Fes 165) Microsoft Office Officel 4 Library: SOLVER\SOLVER XLAM 


Description: Outil d'optimisation et de résolution des equations 


Gérer: |Compléments Excel [=] [atenas 


Cox (amer, 


I.4 Ouvrir Hydrolab2018.xism 


Pour pouvoir utiliser Hydrolab2018 vous devez ouvrir le fichier « Hydrolab2018.xlsm » qui se 
trouve dans le dossier « Hydrolab2018 » directement sous la racine C/ 
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& OS (C:) ® Nom Modifié le 
101_FUJI = 
J 101. i 6) utilisation Hydrolab 201... 21/05/2015 19:06 
à AAA Disque dur JPL = - 
= EF) Hydrolab2010.xis 21/05/2015 18:24 
à AdwCleaner > 
; À G) Hydrolab 2010.docx 13/05/2015 18:44 
d Données pluviométriques d'Algérie , J ; 
Es] Amortissement d'une cr... 13/05/2015 15:02 
à downloads si g 
na | EP) Averse -Cruexls 13/05/2015 15:00 
H Fichiers à traiter = OT ` 
= Ej Q10 méthodes empiriqu... 22/03/2015 21:05 
d) France Inter - 2000 ANS D'HISTOIRE SA sn 
s B) Tab Regmul.xis 18/02/2015 10:01 
d Hydrolab2010 Rues 
E) Tab Anomalie.xls 30/09/2014 19:51 


P" Intel 


SUR — . — 


S’ouvre alors la feuille « Hydrolab2018.xls » qui contient les commandes prêtes à être exécutées 
ainsi que la langue utilisée à l’affichage. Pour changer cette langue ou en ajouter une autre, voir 


plus loin. 
Xi o O T 
Accueil Insertion Mise en page Formules Données Révision Affichage Développeur Compléments 
HYDROLAB2010 ~ 
Commandes de menu 
L1C1 -e _ | Français 
tai Hydrolab2010.xis [Mode de compatibilité] (=) (aj bd 
2 a 
Français JFrançais En= 
2 Ajustement à une loi de Gauss Ajustement à une loi de Gauss 
3 Moyenne= Moyenne= 
4 Ecart-type= Ecart-type= 3 
5 ___Taillen= > pe: Taille n= ___ 
4 4> H| Feuili Feulz < Feu #3 Mafu] » 1 


En activant le menu complément vous devriez voir apparaitre le menu déroulant 
HYDROLAB2018. Vous pouvez alors laisser cette fenêtre ouverte ou la fermer, les macros 
commandes pilotées par le menu déroulant HY DROLAB2018 restent actives. 


H. ELEMENTS COMMUNS AUX STATISTIQUES 


UNIDIMENSIONNELLES 


Sous cette rubrique sont regroupés les ajustements d’une variable aléatoire unidimensionnelle. 
Neuf types de loi ont été retenus : 

Loi exponentielle 

Loi racine-normale 

Loi des extrêmes généralisée (GEV) 


Loi de Gumbel Loi log-normale (Galton) 
Loi normale (Gauss) Loi de Poisson 
Loi de Pearson 3 Loi des fuites 


II.1 Sélection de la plage de données à traiter 


Les données à utiliser sont dans un de vos fichiers Excel. Ce sont aussi bien des valeurs que du 
texte ou des formules. Votre fichier de données sera préservé et les données à traiter seront 
transférées et mises en forme dans un des fichiers « Tab XXX.xlsx ». Le format de vos données 
sera reproduit dans « Tab XXX.xlsx ». (nombre de décimales, formats personnalisé...) 

Les observations de départ doivent être rangées sur une ou deux colonnes, chaque ligne correspond 
à une observation. S’il y a deux colonnes, la première donne le nom des observations (par exemple 
l’année). La première ligne peut éventuellement contenir la nature des observations et des mesures. 
HYDROLAB2018 accepte également que les observations soient rangées sur une ou deux lignes ; 
tout ce qui est dit plus loin reste alors valable en permutant le rôle des lignes et des colonnes 
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IL.1.1 Données en colonne sans nom de variable ni d’observation 


Ici on a sélectionné les 19 hauteurs de pluie annuelles de la station de Contes et on s’apprête à 
ajuster à une loi de Gauss (loi normale). Hydrolab ne traitera que les 18 valeurs numériques et 
1975 sera considéré comme une donnée manquante. 

Une fenêtre s’ouvre alors vous montrant que les données sont en colonne et vous demandant si la 
première ligne contient le nom de la variable. Comme ce n’est pas le cas, on ne coche rien et on 


clique sur OK. 
Lil = 6.» © 3 


| FHER | ACCUEIL INSERTION MISE EN PAGE FORMULES DONNÉES RÉVISION AFFICHAGE DÉVÉLOPPEUR COMPLÉMENTS $ 
Statistiques uni-dimensionnelles + ! + Lo: Exponentelle 
Données multi-dimensonellés + | > Loi de Gumbel 
Commandes de menu 
+ Loilognommale 
F3 z fiine 
+ Loi de Porsson 
A B C + LoiRacine-normale G H i J K 
1 |Cude 6004002 60 + LoiGEy 14001 6048001 6050001 6057001 6071001 6073001 6075001 
2 |STATION ANTIBES 4 CANNE SUR CONTES COURSEGOIL'ESCARENE GUILLAUMEISOLA LEVENS 
3 1967 358.3 > _LoiPesron3 6 7485 6702 6125 7782 7726 
4 1968 1060.6 1162.8 991.9 1634.6 1063.5 1088.1 1177.8 1147.7 
5 1969 768.1 892.5 825.9 1345.1 856 9516 8335 10231 
6 1970 4821 673 8662 PSC SR Ek 
1 1971 8597 956.2 943.8 PR CR ARE RES a 
8 1972 1150.1 13554 11756 G Les valeurs prises por La variable correspondent aux dfférentes lignes 
L} 1973 706.1 725.4 7731 r L 
10 1974 1782 906 611 A 
11 1975 10754 1092 1038.7 
12 1976 1033.7 1156.1 11821 
13 1977 986.5 1031.3 1358.6 
14 1978 8294 7893 6815 
15 1979 1099.4 11736 12925 
16 1980 6333 667 4 584 9 
17 1981 6191 7298 3462 
18 1982 571.5 652 1 709 7 ares | 
19 1983 639.1 693.4 666.9 
20 1984 984 903.5 977 1003.3 1063.4 1029.4 1210.4 1104.9 
pal 1985 6656 6578 5925 975 690 1 729 7106 7296 
2 1986 686.8 7205 737.9 957.2 720.5 687.8 9749 7147 
23 1987 9779 10037 9776 15575 1140.8 967.7 1112.1 1076.3 
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IL.1.2 Données en colonne avec nom de variable et sans nom d’observation 


H ©- e = Données pour ACP.xlsx - Excel 
ACCUEIL INSERTION MISE EN PAGE FORMULES DONNÉES RÉVISION AFFICHAGE DÉVELOPPEUR 


Statistiques uni-dimensionnelles + | > Loi Exponentielle 
Données multi-dimensionelles > | 3 Loi de Gumbel 
Commandes de menu 
> LoilLognormale 
F2 7 f > Loi Normale 
+ Loide Poisson 
A B c > Loi Racine-normale Barre de eree G H | 
1 |Cude 6004002 60 > Loi GEV 14001 __ 6048001 6050001 6057001 6071001 £ 
2 [STATION ANTIBES4 CANNE , CON COURSEGOIL'ESCARENE GUILLAUME“ ISOL 
3 1967 358.3 > LoiPearson3 542.6 748.5 670.2 612.5 
4 1968 1060.6 1162.8 991.9 1260.5 1634.6 1063.5 1088.1 
5 1969 768.1 892.5 825.9 1345.1 856 951.6 
£ Structure de données univariées X 2 Tees es 
8 | |© Les valeurs prises par la variable correspondent aux différentes lignes 2080.7 12133 1161.2 
9 a 1121.8 832.2 815.4 
10 1158.1 745.8 730.2 
11 1800.5 1224.9 955.3 
12 1695 2 1037 1 13121 
13 1747.2 1286.6 1347.5 
14 1158.1 859.7 811.6 
15 1794 1385 1264.2 
16 936.7 794.7 lac lac 
17 1177.9 883.1 lac lac 
18 951 688.7 836.4 
19 | 1070.9 876.4 787.5 
20 1304 304 3U2.9 I 1149.3 1669.3 1083.7 1025.7 
21 1985 665.6 657.8 592.5 707.6 975 690.1 729 


Ici on a sélectionné les 19 données de pluie, mais aussi ne nom de la station « CONTES ». Pour 
que CONTES soit considéré comme le nom de la variable (et non une donnée manquante) il 
faudra cocher la case « la première ligne donne le nom de la variable avant de faire OK. 


IL.1.3 Données en colonne sans nom de variable mais avec noms d’observations 


Cette fonctionnalité nécessite que les données soient dans deux colonnes contigües et que la 
première contienne les noms des observations. Hydrolab2018 a reconnu la présence des noms des 
observations dans la première colonne. Par contre il ne peut savoir si la première ligne contient les 
noms. On laissera ici la première des cases décochée pour taper OK. 
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HS- e-s 
ACCUEIL INSERTION MISE EN PAGE FORMULES DONNÉES RÉVISION AFFIC 


Statistiques uni-dimensionnelles + | + Loi Exponentielle 
Données multi-dimensionelles > | > Loi de Gumbel 
Commandes de menu = 
> Loilognormale 
A3 $ Fie 
> Loide Poisson 
A B c > LoiRacine-normale F | G | 
1 [Code 6004002 602 > LoiGEV 14001 6048001 6050001 Gl 
2 Année ANTIBES 4 CANNE 7 SUR CONTES  COURSEGOIL'ESC 
3 358. 3 > LoiPearson 3 542 6 526.3 7485 
4 _ Structure de données univariées 
6 
7 
: [7 La première colonne donne le nom des observations 
10 
11 
12 
13 


829.4 


15 | 1979 1099.4" 
16 | 1980 633.3 
17 | 1981 619.1 
18 | 1982 5715 


IL.1.4 Données en colonne avec nom de variable et avec noms d’observations 


Cette fonctionnalité nécessite que les données soient dans deux colonnes contigües et que la 
première contienne les noms des observations. On suppose que la première ligne contient le nom 
de la variable. On cochera ici les deux cases avant de taper OK 
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H © = Dor 
ACCUEIL INSERTION MISE EN PAGE FORMULES DONNEES RÉVISION AFFICHAGE DÉVELOPPEUR COMPLÉMENTS 


EEN 


Statistiques uni-dimensionnelles + 


4 


Loi Exponentielle 


Données multi-dimensionelles + Loi de Gumbel 


Commandes de menu 
> Loilognormale 


A2 X Š f J Loi Normale 


> 

> Loide Poisson 
> LoiRacine-normale F G H | J 
> 


A B c K 

1 [Code 6004002 60! > LoiGEV D 273001 607500 
2 5 ANTIBES À ACANNE x LEVENS 

3 CT > LoiPearson3 778.2 7721 

4 1162.8 991.9 177.8 1147. 

5 892.5 825.9 333.5 1023. 

6 673 868.2 1137 1152: 

7 956.2 943.8 390.4 1390. 

8 1355.4 1175.6 123.5 1534. 

9 725.4 773.1 791.7 880.: 

906 611 597.3 759- 

1092 1038.7 235.8 12221 

1156.1 1182.1 270.1 11631 

1031.3 1358.6 446.4 1357. 

789.3 681.5 301.8 982. 

1173.6 1292.5 322.4 1452: 

667.4 584.9 673.7 676.9 936.7 794.7 lac lac 795.1 


IL.1.5 Lancement du traitement 


En cliquant sur le bouton OK, vous lancer le traitement statistique souhaité. Vos données restent 
préservées dans le fichier d’origine mais les valeurs en sont copiées dans un fichier Tab XXX.xlsx 
qui s’ouvre directement pour faire apparaître les données et les résultats de leur traitement. Ces 
fichiers sont protégés en écriture. Si vous souhaitez néanmoins conserver une feuille de résultats 
enregistrez-la sous un nouveau nom. 


IL2 Les résultats numériques 


Après avoir lancé un ajustement vous verrez s’ouvrir une nouvelle feuille contenant vos données 
et les résultats sous forme numérique et graphique Vous pouvez y changer les formats, les 
polices... Si vos données sélectionnées avaient un format particulier, il sera reproduit dans cette 
feuille. Toutes les feuilles (à l’exception de la loi de Poisson) se présentent de façons analogues et 
vous y trouverez : 


justement à une loi de Gauss 
l Moyenne= 835.4818182 
-ty 196.2031744 

V z Jari 


classement 

1 386.2840077 | 587.2137980 
2 $ 620.2854776 | 527.6643657 | 690.8130409 
3 0.2273 .7 6 688.7495879 | 608.2990138 | 754.1356137 
4 0.3182 47 709. 742.7190918 | 669.9902709 | 805.9242479 
5 0.4091 è 790.3778246 | 722.7562169 | 853.3687434 
6 0.5000 . 766. 835.4818182 | 771.0216295 | 899.9420068 
7 0.5909 ? À 880.5858118 | 817.5948929 | 948.2074194 
8 0.6818 4 b 928.2445445 | 865.0393884 | 1000.973365 
9 0.7727 .7 i 982.2140485 | 916.8280227 | 1062.664623 
10 0.8636 .097 ë 1050.678159 | 980.1505955 | 1143.299271 
11 0.9545 j 1167.187149 | 1083.749838 | 1284.679629 
Lac. 


11.2.7 


IL2.1 Les paramètres ajustés de la loi de distribution 


Ici, pour la loi de Gauss, il n’y a que deux paramètres : la moyenne et l’écart-type. 


IL.2.2 La taille de l’échantillon sélectionné et le nombre de données utilisables 


Dans ce cas de figure on avait sélectionné douze observations mais une était lacunaire 
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IL2.3 Le test d’Anderson 
Pour chaque ajustement on calcule la valeur du Wn2 d’Anderson-Darling : 


wW? = n-1S (2i 1)Ln[F(x, )]+(2n-2i+1)Enll-F(x;)|]) 
ni 


Et on évalue la probabilité pour que cette valeur soit dépassée du seul fait du hasard alors que la 
loi ajustée serait la bonne. Ici la probabilité d’être dépassée est de 0.295 (29.5%) cette probabilité 
est suffisamment forte pour ne pas rejeter l’hypothèse d’un ajustement satisfaisant. A priori on 
retient les seuils suivants : 

- proba >0.1 (10%) ajustement satisfaisant (surligné de vert) 

-  0.1(10%)> proba >0.05 (5%) ajustement douteux (surligné d’orange) 

- proba <0.05 (5%) ajustement à rejeter (surligné de rouge) 


IL.2.4 Les intervalles de confiance 


A priori on retient les intervalles de confiances à 70% , c'est-à-dire que pour chaque fréquence 
théorique on évalue les bornes inférieures et supérieures entre lesquelles il y a 70% de chance que 
se trouve le véritable quantile. Il est tout à fait loisible de changer cet intervalle de confiance par 
exemple en tapant 90 au lieu de 70. Tout se remet à jour dans les résultats numériques, les 
graphiques et les tableaux d’application probabiliste. 


IL.2.5 Les noms des observations classées 


On trouve dans cette colonne les noms des observations classées dans l’ordre des valeurs 
numériques observées. Ici la plus petite observation est la première et la plus grande, la onzième. 
La 9% observation est lacunaire. Si on avait également sélectionné la colonne contenant les années 
d’observation on lirait directement que la plus petite observation est celle de 1967 et la plus grande, 
celle de 1977. L’année 1975 est lacunaire. 


IL2.6 Les valeurs numériques observées classées par ordre croissant 


IL2.7 Les fréquences expérimentales et les variables réduites associées 


Pour chaque valeur classée xi, on évalue sa fréquence expérimentale par : Fexp(xi)= (1-0.5)/n 
A chaque fréquence expérimentale est associée une valeur de la variable réduite adaptée à la loi. 


IL2.8 Valeurs expérimentales, théoriques et bornes de l’intervalle de confiance 


Pour chaque variable réduite (et donc chaque fréquence expérimentale) on donne la valeur 
observée, la valeur théorique selon l’échantillon et les bornes entre lesquelles il y 70 % de chance 
que se trouve la valeur observée. 


IL3 Le graphique d’ajustement 


Dans la partie droite de la feuille vous trouverez un graphique d’ajustement que vous pouvez 
modifier selon votre gout et sur lequel vous pouvez ajouter des commentaires libres. Ce graphique 
est un objet que vous pouvez copier et coller dans n’importe quel document word ou autre. 


IL3.1 En abscisse la variable réduite associée aux fréquences 
Cette échelle dépend de la loi choisie. 


IL3.2 En ordonnée les valeurs de la variable aléatoire ou de sa transformée (log, racine...) 


En général les abscisses et ordonnées sont choisies de telle façon que la loi théorique soit 
représentée par une droite. Vous pouvez librement choisir les graduations et leurs limites 
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Ajustement à une loi de Gauss 


1400 - 11.3.4 (Moyenne=835 48, Ecart-type=196.20. Taille n=11, LC. à (en%)=70) 


1200 - 


1000 - 


11.3.1 


£ 

"2800 

Z 

Li 

g 

#600 oints expérimentaux 

G . z . . z 

> loi théorique ajustée 11.3.3 


intervalle de confiance 


200 - 


-2.5 2 -15 -0.5 0.5 1.5 2 


Fréquences selon une échelle gaussienne 


IL3.3 Courbe des valeurs expérimentales, théoriques et bornes de l’intervalle de confiance 


Vous y trouverez en nuages de points les valeurs expérimentales et en courbe grasse les valeurs 
théoriques et en courbes fines, les bornes de l’intervalle de confiance. 


IL3.4 Récapitulatif de l’ajustement 


Dans le graphique sont figurés les éléments numériques essentiels de l’ajustement. En général les 
abscisses et ordonnées sont choisies de telle façon que la loi théorique soit représentée par une 
droite. Vous pouvez librement choisir les graduations et leurs limites 


IL4 Les tableaux d’exploitation probabiliste 


A la fin de la feuille de résultats vous trouverez deux tableaux tels que les suivants : 


Variable Valeur Borne Borne Période de 
réduite théorique inférieure supérieure Valeur Fréquence retour 


Celui de gauche permet d’évaluer la valeur de la variable aléatoire et son intervalle de confiance 
pour une fréquence donnée. Celui de droite permet à l’inverse, d’associer à une valeur particulière 
de la variable aléatoire, sa fréquence et sa période de retour. 

A titre d’exemple il suffit de taper à la place du point d’interrogation la valeur 0.9 pour savoir qu’à 
cette fréquence la valeur de la variable réduite est de 1.282 (ici variable réduite de Gauss) 


On constate que pour cette fréquence Variable Vs Bass Biss 
Fréquence réduite théorique inférieure supérieure 


la valeur xr la plus probable de la 
| o9 | 1282 | 1086926303 | 1012.843913 | 1186.823728 


variable aléatoire est de 1086.9 et 
qu’il y a 70% de chance d’avoir : 
1012.8< xr <1186.8 
De la même façon, on peut évaluer pour n’importe quelle valeur, sa fréquence théorique et sa 
période de retour en remplaçant le point d’interrogation par la valeur désirée. 
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Ainsi une pluie de 1126 mm correspond à une 
fréquence théorique F de 0.931 soit une période de 
retour T=1/(1-F)-14.4 années. Attention ceci n’a de 
sens que si l’on travaille sur un échantillon constitué 
de valeurs annuelles. 

Si vous désirez faire ces calculs pour plusieurs valeurs il suffit de recopier la dernière ligne des 


tableaux vers le bas comme l’illustre le schéma ci-dessous. 
Variable Valeur Borne Borne Période de 
réduite théorique inférieure supérieure Valeur e retour 


1086.926303 | 1012.843913 | 1186.823728 


1158.207321 | 1075.900085 
1238.433874 | 1145.461351 | 1372.776249 


ME E E 


HI. SPECIFICITES DES DIFFERENTES LOIS 


HI.1 Loi de Gauss (ou aussi loi normale) 


IL.1.1 Domaines d’applications de la loi de Gauss 


Cette loi est très fréquemment utilisée en hydrologie. Son emploi se justifie principalement par le 
théorème central limite qui établit que la somme x de p réalisations indépendantes d’une variable 
aléatoire quelconque y tend vers une loi de Gauss lorsque p tend vers l’infini. Ainsi cette loi 
s’applique assez bien à des variables hydrologiques « moyennes » telles que la pluviométrie 
annuelle, les débits moyens annuels... 


IIL.1.2 Compléments techniques 


La mise en œuvre de la loi de Gauss correspond intégralement au dispositif commun de lancement 
et de présentation des résultats. Cette loi est définie sur [-c,++|et peut convenir à n'importe quel 
échantillon. Il nous +. seulement bon de préciser quelques points techniques. La fonction de 


répartition F(x) = e 2duavec u=% X dépends de deux paramètres qui sont la moyenne x 
(o) 


et l'écart-type o. 
e Ces paramètres moyenne et écart type sont évalués selon les estimateurs sans biais suivant : 


i=1 
n 
e L'intégrale de Gauss est évalué par les fonctions LOINORMALE() et 
LOI.NORMALE.INVERSE() 


X = 


, fonctions MOYENNE() et ECARTYPE() d’Excel 


e pour l'intervalle de confiance à a % sur le quantile xF, les bornes sont évaluées par : 
XF + Erreur !6 (Informations Techniques du CTGREF, Cahier 31, N°2, 1978) 
(2 fois plus pour la borne supérieure et 2 fois moins pour la borne inférieure) 
tF : variable réduite de Gauss ayant la fréquence au non-dépassement F 
; a r 7 l-a 
tą : variable réduite de Gauss ayant la fréquence au non-dépassement 1- =- 


e pour l'intervalle de confiance à a % sur les moyennes et écart-type (Tableau au-dessus ou sous 
le graphique) on utilise classiquement les fonctions LOI.STUDENT.INVERSE.N() et 
KHIDEUX.INVERSE(). 
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II.2 Loi de Galton (ou aussi loi log-normale) 


IIL.2.1 Domaines d’applications de la loi de Galton 


Dans Hydrolab la loi de Galton est utilisée sous une forme simplifiée où l’on suppose que c’est le 
logarithme de x qui suit une loi de Gauss. Son emploi se justifie pour des grandeurs hydrologiques 
résultants de la multiplication de nombreux facteurs. En effet si les facteurs y se multiplient les 


log(y) s’ajoutent : x= lI y, => log(x) = D log(y,) et log(x) tends ver une loi de Gausse 
i=1 i=1 


lorsque n augmente. 

Ceci explique que cette loi soit souvent utilisée dans la modélisation de phénomènes extrêmes tels 
que crues ou étiages. 

Dans certains cas de figure, la variable x n’est pas bornée par zéro. Des prélèvements en eau, des 
infiltrations ... peuvent justifier l’existence d’une borne inférieure Xo à la variation de x. C’est 
alors le log(x-xo) qui suit une loi de Gauss et on supposera que ce seuil est parfaitement connu et 
non ajusté. 


II.2.2 Compléments techniques 


L’ajustement à une loi de Galton reprends toutes les étapes de l’ajustement à une loi de Gauss mais 
en travaillant sur le log(x) puis en revenant aux valeurs naturelles de x. Si l’on se donne une borne 
inférieure xo, il suffit de la rentrer manuellement et on travaille alors sur le log(x-xo) puis on revient 
aux valeurs naturelles de x. Mais attention, les intervalles de confiances supposent alors que xo 
est parfaitement connu. 


IIL3 Loi des fuites 


IIL3.1 Domaines d’applications de la loi des fuites 


La loi des fuites préjuge que la grandeur x est la somme de k réalisation d’une grandeur y. On 
suppose que le nombre k de réalisation suit une loi de poisson, et que la grandeur y suit une loi 
exponentielle de paramètre de position nul. Ces hypothèses sont cohérentes avec de nombreuses 
variables hydrologiques telles que les précipitations mensuelles et même annuelles. 


IIL.3.2 Compléments techniques 


Dans le logiciel Hydrolab, nous nous sommes inspiré essentiellement des travaux de Pierre 
RIBSTEIN (Cahier ORSTOM série hydrologie Vol 20 N° 2, 1983) qui s’inspire lui-même de la 
thèse de C. BARBUSIAUX (Etude statistique de la loi des fuites, thèse de 3eme cycle, Facultés 

des sciences de Paris, 1969). 
e Sile nombre k de réalisation de y suit une loi de Poisson de paramètre de forme 1 : 

Ex 

k! 

e Si la grandeur y suit une loi exponentielle simple de paramètre d’échelle s: 


Fo)=1-e% 


Proba (k) = 


k 
La grandeur x = >. y, suit une loi des fuites de fréquence au non dépassement Fx) : 
i=0 


PE ¿į u’ x 
3 — avec u=— 
S 


F(x)=1-e% e* - - 
m (+1) J! 
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œ 50 

Dans la feuille de calcul et pour le calcul de F(x) la somme J, est limitée à D , ce qui est 
i=0 i=0 

largement suffisant pour que les termes manquants soient bien négligeables. 

Pour l’inversion de cette fonction, nous avons calculé des couples x — F(x) en nombre suffisant 

puis interpolé entre les points. 

Pour l’ajustement des paramètres de forme 4 et d’échelle s, nous avons retenu la méthode 

proposée par C. BARBUSIAUX. Elle donne des estimations fiables de ces paramètres en se basant 

sur le nombre no de valeurs nulles dans l’échantillon de taille n, sur la moyenne y et la variance 

V: 


X xX 


E a ) 


FO If) Li) 
n n 


2: 
200860) Ca) te 
n yV n n 

Il n’existe pas d’étude exhaustive des intervalles de confiance sur les estimations des paramètres 
et des quantiles. Nous avons donc procédé aux tirages de 1000 échantillons de taille 12, 25 50 et 
100 pour des valeurs de À de 0.5, 0.75, 1, 1.5, 2, 3, 4.5, 7 et 10. Nous en avons extraits les 
intervalles de confiance à 50%, 70%, 80%, 90% et 95% sur les paramètres 2 et s et sur les 
quantiles de Fréquence 0.001, 0.002, 0.005, 0.01, 0.02, 0.05, 0.1, 0.2, 0. 5, 0. 8, 0. 9, 0. 95, 0. 98, 
0.99, 0.995, 0. 998 et 0. 999 ! C’est à partir de cette information qu’Hydrolab interpole les 
intervalles de confiance. 

Pour les paramètres 2 et s, l'intervalle de confiance est donné dans le tableau au-dessus ou sous 
le graphique d’ajustement. 


II.4 Loi racine normale 


IIL.4.1 Domaines d’applications de la loi racine normale 


De façon un peu analogue à ce qui a été dit pour la loi de Galton, on suppose ici, que c’est la racine 
carrée de x qui suit une loi de Gauss. Il n’existe pas de justification très pertinente pour l’utilisation 
de cette loi. Cependant elle est d’un emploi fréquent car elle est une approximation assez correcte 
et simple à traiter, de la loi de distribution des fuites que nous venons de voir 

On constate que certaines grandeurs hydrologiques sont la somme de réalisations indépendantes, 
comme nous l’avons dit pour la loi de Gauss. Cependant si leur nombre ne tends pas vers l’infini, 
cette distribution présente une dissymétrie positive. C’est le cas pour des précipitations mensuelles 
mais aussi annuelles si la pluviométrie est faible et avec peu d’événements pluvieux dans l’année. 
On peut également remarquer que même si les pluies annuelles sont fortes en moyenne, mais avec 
un écart-type faible, le passage de x à 4/x est quasi linéaire. La loi racine normale parait donc très 
bien adaptée pour toutes les pluies annuelles sur le bassin méditerranéen. 


IIL.4.2 Compléments techniques 


L’ajustement à une loi de racine normale reprends toutes les étapes de l’ajustement à une loi de 
Gauss mais en travaillant sur4/x puis en revenant aux valeurs naturelles de x. On relèvera 
néanmoins la mauvaise adéquation de cette loi, si les valeurs de x sont fréquemment nulles. En 
effet la loi de Gauss est définie sur [-c,+«] alors que ~x est positif ou nul. Avec Hydrolab les 
valeurs négatives de Vx sont ramenées à zéro mais ce n’est qu’un pis-aller ! 
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IILS Loi de Gumbel 


HI.5.1 Domaines d’applications de la loi de Gumbel 


Cette loi a longtemps été considéré comme « la » loi de distribution des extrêmes. Comme on le 
verra plus loin, ce n’est en fait, qu’une des trois formes possibles de distributions des extrêmes. 
Cependant cette loi présente l’énorme avantage de ne dépendre que de deux paramètres et 
également de disposer d’estimateur simples des paramètres. Ceci explique que cette loi soit très 
fréquemment utilisée pour l’étude des crues maximales annuelles, des précipitations journalières 
maximales annuelles … 


II.5.2 Compléments techniques 


X— Xo 


La fonction de répartition est : F(x)=e°" avec u= (variable réduite de Gumbel) On notera 


la simplicité du passage des valeurs de x à leurs fréquences et réciproquement. 

Le terme u est la variable réduite de Gumbel ; Xo est le paramètre de position (mode) et g est le 
paramètre d'échelle différent de zéro et positif (g est aussi appelé "gradex"). 

Cette loi est définie sur |- œ,+%0]et peut convenir à n'importe quel échantillon. 

On peut noter dès à présent le comportement asymptotiquement exponentiel de la distribution de 
Gumbel : si F tend vers 1, en posant T = Erreur !, la variable réduite de Gumbel u tend vers Ln 
T. 


Les paramètres sont estimés par la méthode des moments : g = 0,78 oet x,= xi - 0,577 s 
L'intervalle de confiance à a % sur un quantile xp s'exprime en fonction de l'écart-type © par : 

?r -h]O<XF< ik +h2 6 
où h1 et h2 sont des paramètres dépendant de la taille n de l'échantillon de la fréquence F et de la 


valeur de a. 
e hı et h2 seront évalués par la formule suivante (avec le signe + pour h2 et le signe - pour h1) : 


2 
u u 
a +113, +112 E E E) 
h > = n 7 (Informations Techniques du CTGREF, Cahier 31, N°2, 1978) 
u 
1-11 
n 


; p Š ; ; 1- 
e u, est la variable réduite de Gauss correspondant à la fréquence au non-dépassement Eee 


e tp est la variable réduite de Gumbel correspondant à la fréquence au non-dépassement F, ramenée 
à sa moyenne et à son écart-type : tp = Erreur ! 


IIL6 Loi exponentielle 


IIL.6.1 Domaine d’applications de la loi exponentielle 


La loi exponentielle s’applique assez bien aux observations de valeurs supérieures à un seuil. C’est 
souvent le cas pour les pluie journalières supérieures à un seuil, les débits de crue supérieures à un 
seuil...Elle se justifie aussi indirectement en remarquant que si les valeurs supérieures à un seuil 
suivent une loi exponentielle, et si le nombre de dépassement du seuil suit une loi de Poisson, alors 
les valeurs maximales annuelles suivent une loi de Gumbel. 


IIL.6.2 Compléments techniques 


La fonction de répartition est : F(x) =1 Le £ 


Utilisation d'Hydrolab 2018 -20- 
POLYTECH 
NICE-SOPHIA Décembre 2017 


Cette loi a donc deux paramètres (x, et g) et est définie sur l'intervallefx,,+%]. Nous considérons 
ici que x, est la borne inférieure de l'intervalle de définition et qu'il est donc connu a priori. Dans 
la feuille la valeur par défaut de x, est la plus petite valeur de l'échantillon mais on peut toujours 
se donner une autre valeur de x, en la tapant directement dans la cellule située à droite de " x=". Il 
faut bien sûr, que x, soit inférieur ou égal à la plus petite valeur de l'échantillon. 

Le paramètre g (que l'on peut également appeler gradex) est estimé par :g=x-x, 

L'intervalle de confiance à a % sur un quantile xp s'exprime ainsi : 


ka, Pen) «x, <ôp tuy Pet 


; 7. i : z J= 
e u, est la variable réduite de Gauss correspondant à la fréquence au non-dépassement Eee 


(formule simplifiée d'après MIQUEL 1984 et valable si g est le seul paramètre ajusté). 


II.7 Loi des extrêmes généralisée (GEV ou loi de Jenkinson) 


IIL.7.1 Domaine d’applications de la GEV 


Jenkison (1955) a montré que les lois de distribution des extrêmes pouvaient se mettre sous une 


forme unique : 
Jeo] 
S 
F(x)=e 


Trois paramètres interviennent : Xo le paramètre de position, s le paramètre d’échelle et k le 
paramètre de forme. Sous cette formulation unique on retrouve en fait les trois types de distribution 
selon les valeurs prise par k : 

Si k est positif, on retrouve une loi proposée par Weibull (GEV de type II) où x est défini sur 
l’intervalle ]-,x,'[ Cette borne supérieure de l’intervalle de définition se déduit des paramètres 
s, k et Xo par la relation : Xo=p+X0. En fait ce type de distribution se rencontre rarement en 
hydrologie et on la réserve à des études de valeurs extrêmes minimales (température, étiages...). 


LATE Gen) 
Si k est égal à zéro on constate que j-e >e s 


. On retrouve donc la loi de Gumbel 
(GEV de type I) vue précédemment et x est défini sur l’intervalle ] -œ ,+oo [ 


Enfin, si k est négatif, on trouve une loi proposée par Frechet (GEV de type ID) où x est défini sur 
l’intervalle ]x}',+œ[ Cette borne inférieure de l’intervalle de définition se déduit des paramètres 


s, k et xo par la relation : X0 = +X0 Ce type de loi se rencontre souvent en hydrologie (cues, 


précipitations extrêmes.…..), mais la difficulté va résider dans l’estimation du paramètre de forme 
qui est très soumise aux aléas de l’échantillonnage. Il parait très hasardeux de caler la valeur de k 
sur un seul échantillon. Il est plus raisonnable de choisir a priori la valeur de k d’après une étude 
régionale. Koutsoyiannis (2004) a étudié 169 séries longues de pluie extrêmes (Etats-Unis, 
Royaume-Uni, France, Italie, Grèce) et il propose de fixer k à -0.15. Pour les pluies extrêmes 
d'Algérie (508 séries de 43 années en moyenne), je propose de retenir k=-0.073. 


IIL.7.2 Compléments techniques 


Dans le cas général, il est préconisé d’utiliser la méthode des moments pondéré (Voir LUBES & 
MASSON, Hydrologie Continentale, vol.6 n°1, 1991). On estimera tout d’abord les trois premiers 
moments pondérés bo, bı et b2 : 
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hl b=15 i—] -15 GDG2) , 
D > Le Dire Iaa 


(n étant la taille de l ehant Ton eti : rand dans l'échantillon classé par ordre croissant) 
On peut montrer que ces trois premiers moments pondérés bo, bı et b2 , sont liés au trois paramètres 
Xo, S et k par les relations : 

3b,-b, 1-3 * 
2b, -bọ 1-2™ 


, bo =xo+pli-TU+K)] et 2b; bo =ETU+K a+ 


La fonction gamma est telle que : r(1 +k) = [ue "au et on peut aisément montré que 


Pour obtenir la valeur de F (k + 1) sur Excel, il suffit de faire appel à la fonction GAMMA(). 

La première équation ne contient qu’une seule inconnue, le terme k, et on peut la résoudre par 
itérations successives. Dans le cas général où —-0.5<k<0.5 on peut évaluer k explicitement par la 
2b,-b,  Ln(2) 

3b,-b, aG) 

On en déduit immédiatement les valeurs de s et x, par : 


2b, -b,)k T(1+k)-1 
5 ne e Xo =b, 3 POEREX 
{-2*)T(+k) k 
Il faut être très prudent sur cette estimation de k et nous recommandons de choisir une valeur 
régionale en la tapant directement dans la case concernée. 


Pour évaluer les intervalles de confiance sur l’estimation des quantiles, nous nous sommes 


relation : k=78590c+2.9554c2 avec c= 


contenté d’envisager uniquement les lois de Fréchet. En posant Xo'= +0, la variable 


y= Ln(x-x,) suit une loi de Gumbel. Nous avons donc calculé les intervalles de confiance sur 
les quantiles de y puis nous sommes revenus aux valeurs de x. Ceci n’est qu’une approximation 
qui suppose que l’on ne se trompe pas sur la borne x; 


II.8 Loi de Pearson III ou loi Gamma incomplète 


IILS8.1 Domaine d’applications de la loi de Pearson III 


Il n’y a pas de justification théorique à l’utilisation de cette loi en hydrologie, cependant elle est 
utilisée par pragmatisme dans certains secteurs. Par exemple on constate que les débits moyens 
annuels d’Algérie s’ajustent assez bien à ce type de loi. 


IIL.S8.2 Compléments techniques 


La loi Pearson III est une loi à 3 paramètres dont la fonction de répartition est donnée par : 


L Popa > X—X 
F(x)=——|u""e"du avec u= 4—2 
FC) 5 s 
L'intervalle de variation de x est [ xo, œ [ et les trois paramètres d'ajustement sont : 
Xo : paramètre de position (borne inférieure) 


s : paramètre d'échelle (de même dimension que x) 
y : paramètre de forme (positif différent de zéro) 
En hydrologie, la borne inférieure xQ est très souvent fixée à zéro. Les deux paramètres restants 
2 2 
s’ajustent alors par la méthode des moments par les expressions s= 2 et A= où o est 
X X 
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l’écart type et x la moyenne de la série. On remarque alors que le paramètre de forme est 
1 

cy? 

L’évaluation de F(x) et son inversion se font avec les fonctions d’Excel LOIGAMMA.NO et 

LOI.GAMMA.INVERSE.N() 

Il n’existe pas d’étude exhaustive des intervalles de confiance sur les estimations des paramètres 

et des quantiles. Nous avons donc procédé aux tirages de 1000 échantillons de taille 10, 20, 50, 

100 et 200 pour des valeurs de CV de 0.2 à 2. Nous en avons extraits les intervalles de confiance 

de 50%, à 99% sur les quantiles de fréquence 0.01, 0.02, 0.05, 0.1, 0.2, 0.3, 0.4, 0. 5, 0.6, 0.7, O. 

8, 0. 9, 0. 95, 0. 98 et 0.99 C’est à partir de cette information qu’Hydrolab interpole les intervalles 

de confiance. 


II.9 Loi de Poisson ou loi des événements rares 


IIL.9.1 Domaine d’applications de la loi de Poisson 


Soit p la probabilité pour qu'un événement se produise lors d'une expérience et q la probabilité 
pour qu'il ne se produise pas ; on a bien évidemment : p + q = 1. Supposons maintenant que l'on 
effectue n expériences indépendantes ; on cherche alors quelle est la probabilité P (k) pour que 
l'événement se produise k fois dans les n expériences. On a alors : 
n! Eak 
Pa © 

Cette expression à deux paramètres (n et p) est connue sous le nom de loi de Bernoulli ou encore 
de loi binomiale. Il est peu fréquent d’utiliser cette loi en hydrologie, par contre on est souvent 
amené à travailler sur les nombres d’occurrence d’évènement rares. Si l’on s’intéresse au nombre 
k de crues dont le débit Q a dépassé un certain seuil Qo lors d’une année on se rapproche du schéma 
de Bernoulli. En effet on pourrait admettre que dans l’année il y a n=365 jours, et que la probabilité 
p pour qu’une crue y dépasse Qo est une valeur très faible voisine de zéro. Donc si n est très grand 


e™” ut 
k! 

encore loi des évènement rares. Elle ne dépend que du seul paramètre y représentant le nombre 

moyen de réalisation dans l’année. 

Le nombre d’événement dépassant annuellement, un seuil comme les crues, les pluies 

journalières... suit bien souvent une loi de Poisson. 


et si q=1-p tend vers un, la probabilité de k devient: P(k)= C’est la loi de Poisson ou 


IIL.9.2 Mise en œuvre 


C’est la seule loi discrète développée dans Hydrolab et sa mise en œuvre est singulière ! Ceci 
mérite de traiter un exemple explicité à la page suivante. 

On a relevé sur 36 années les hauteurs de pluies journalières et leur date d’occurrence. On a ensuite 
classé ces hauteurs et dates par hauteurs décroissantes. On a enfin transformé ces dates en année 
hydrologique débutant au 1° septembre. Par exemple le 13/10/1973 appartient à l’année 
hydrologique 1973-1974 et le 19/3/1979 appartient à l’année hydrologique 1978-1979. Les 42 
années d’apparition des pluies supérieures à 60 mm sont alors sélectionnées (Cellules M6 à M47) 
et on lance dans la rubrique Complément puis Hydrolab2018 puis Statistiques unidimensionnelles 
et enfin loi de Poisson. 

On obtient la feuille Tab Poisson où il suffit de préciser le nombre d’année d’observation (36 dans 
la cellule G3). 
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i | A | B c An] E | F x y l z 

3 lle 5 l Ajustement à une loi de Poisson 

1 FA Nombre de dates rencontrées 

| = = 3 Nombre de dates possibles 

5 Année Date Pj 4| Moyenne 1,17 at-type théo. = 1,08 

6 [1973-1974 | 13/10/1973 191,4 5 

7 [1987-1988 05/10/1987 117,8 G 

8 30/09/1998 116.6 T Dates rencontréeates rencontrée Effectif 

9 24/10/1999 108,2 8 0 13 1121 

10 03/10/1992 104,2 9 1 12 13,08 

11 24/12/2000 98,4 2 7 7,63 

12 10/10/1987 95,8 3 2 2,97 

13 03/11/1968 93,1 4 1 0,87 

14 08/10/1998 90,2 5 0 0,20 

15 |: 22/09/1996 89,8 6 1 0,04 

16 06/10/1992 89 

17 19/12/1997 89 Ajustement à une loi de Poisson 

18 07/10/1993 87,4 

19 19/03/1979 85 

20 29/09/1991 84,8 

21 20/06/1966 80,6 20 

22 19/09/1999 77.4 21 

23 [2000-2001 05/11/2000 76,8 2 

24 14/09/1968 76,3 23 

25 (2000-2001 | 23/11/2000 76,2 24 

26 16/10/1990 75.6 = 

27 28/09/1974 753 z 

28 10/01/1994 73,4 

29 08/12/1990 732 PA FRERES E (Moyenne=1,167) 
29 1966-1967 3 

30 05/10/1993 70,2 r 1968-1969 4 - 

31 15/10/2000 68.6 on + 

32 03/11/1966 66,6 1971-1972 1 

33 25/02/1978 66,1 1973-1974 1 

34 | 1970-1971 13/11/1970 65,7 1974-1975 1 

35 [1966-1967 | 29/09/1966 65,1 19771078 5 

36 1979-1980 | 12/10/1979 63,2 1978-1979 1 

37 (1984-1985 | 03/10/1984 62,8 1979-1980 2 

38 |2000-2001 11/10/2000 62,8 1983-1984 1 

39 {2000-2001 13/11/2000 62,2 J 1984-1985 1 

40 |1988-1989 01/12/1988 62 40 1987-1988 2 

41 |1979-1980 20/09/1979 61,8 41 1988-1989 1 

42 | 1966-1967 15/10/1966 61,4 42 1990-1991 2 

43 | 1968-1969 13/03/1969 61,4 43 1991-1992 1 

44 | 1968-1969 16/11/1968 61,3 44 1992-1993 2 

45 | 1971-1972 19/02/1972 61,1 45 1993-1994 3 

46 [1983-1984 | 22/05/1984 60.6 46 1996-1997 1 

4T (A977-4978 | 12/01/1978 60,2 47 1997-1998 1 

48 | 1976-1977 20/02/1977 60 48 1998-1999 2 

49 |1967-1968 23/08/1968 59 49 1999-2000 2 

50 1982-1983 18/06/1983 59 50 | 2000-2001 6 


IIL.9.3 Compléments techniques 


Dans la feuille de résultats Tab Poisson, on retrouve en colonne A l’échantillon des années 
hydrologiques, dans les colonnes B et C on retrouve pour chaque année où il y a eu au moins un 
dépassement, le nombre de dépassement. Par exemple il y a eu 1 pluie supérieure à 60mm en 1965- 
1966, 3 pluies en 1966-1967 et 4 pluies en 1968-1969... 

Le seul paramètre de la loi de Poisson est le nombre moyen de dépassement (avec 42 pluies en 46 


ans, On à LH = z =1.17 indiqué en cellule B4). Si la distribution du nombre d’occurrence suit bien 


une loi de Poisson, on démontre que la moyenne est lié à l’écart-type ø par la relation : o = Ju 


Donc on devrait avoir un écart-type théorique de 1.08 = V1.17 On constate que l’écart-type 
expérimental est un peu plus fort avec o=132 Sur le graphique joint on a reporté les nombres 
théoriques et expérimentaux des différentes occurrences. Ce graphique permet de se faire une idée 
subjective de la qualité de l’ajustement. 


POLYTECH 
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